Rapid-purification protocols for optical homodyning 
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We present a number of rapid-purification feedback protocols for optical homodyne detection of 
a single optical qubit. We derive first a protocol that speeds up the rate of increase of the average 
purity of the system, and find that like the equivalent protocol for a non-disspative measurement, 
this generates a deterministic evolution for the purity in the limit of strong feedback. We also 
consider two analogues of the Wiseman- Ralph rapid-purification protocol in this setting, and show 
that like that protocol they speed up the average time taken to reach a fixed level of purity. We 
also examine how the performance of these algorithms changes with detection efficiency, being an 
important practical consideration. 
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I. INTRODUCTION 

Rapid-purification protocols increase the rate at which 
the state of a system ispurified by a continuous measure- 
ment 0, S, S i, 1, S UMM- The y do this by applying 
feedback control to the system as the measurement pro- 
ceeds. All such protocols described to date have been 
devised for continuous measurements of an observable 
(that is, measurements that are not dissipative). Un- 
der this kind of measurement the evolution of the system 
density matrix, p, is given by the stochastic master equa- 
tion GaBlEa, 



dp = -(i/h)[H,p]dt-k[X,[X,p]]dt 
+V2k(Xp + pX - 2(X)p)dW, 



(1) 



where X is the hermitian operator corresponding to the 
observable being measured, H is the Hamiltonian of the 
system, dW is Gaussian white noise satisfying the Ro cal- 
culus relation dW 2 = dt. The observers continuous mea- 
surement record, which we will denote by r(t), is given 
by dr = (X)dt + dW/V&k. This kind of measurement 
will project the system onto an eigenstate of X after a 
time t 3> l/(A 2 fc), where A is the difference between the 
two eigenvalues of X that are nearest each other. 

Photon counting and optical homodyning do not fall 
into the above class of measurements because they sub- 
ject the system to dissipation. Thus if one has a single 
optical qubit, consisting of a single mode containing no 
more than one photon, and one measures it with a pho- 
ton counter, then regardless of whether the measurement 
tells us that the state was initially |0) or |1), as t — > oo 
the final state is always |0). If we wish we can think of 
this as a measurement of the photon number (that is, a 
measurement in the class above with X = a'a), followed 
by an irreversible operation that takes both |0) and |1) 
to the vaccum. 

Our purpose here is to examine whether there exist 
rapid-purification feedback protocols for homodyne de- 
tection performed on a single optical qubit, and if so, to 
compare their properties with those pertaining to a con- 



tinuous measurement of an observable on a single qubit. 
Our motivation is partly theoretical interest regarding 
the effect of dissipation on rapid-purificaton protocols, 
and partly to explore whether such protocols can be im- 
plemented in an optical setting. Before we begin it is 
worth recalling the properties of the single-qubit rapid- 
purification protocols that have been derived to date for 
non-dissipative measurements. The first is the protocol 
introduced by one of us (see also ||) in which one 
applies feedback control to speed up the increase in the 
average purity of the system. The average here is taken 
over all possible realizations of the measurement (all pos- 
sible measurement records r(t)). The protocol involves 
applying feedback during the measurement to keep the 
Bloch vector of the state of the qubit perpendicular to 
the basis of the measured observable, X. In the limit of 
strong feedback, and high final average purity, this pro- 
vides a factor of two decrease in the time required to 
reach a given average purity. In the limit of strong feed- 
back the protocol also eliminates the stochasticity in the 
purification process, so that the purity increases deter- 
ministically. 

The second protocol, introduced by Wiseman and 
Ralph (see also 0, Q), involves applying feedback 
to keep the Bloch vector parallel to the basis of the 
measured observable. (If the system has no apprecia- 
ble Hamiltonian, then the measurement will do this of 
its own accord, and feedback is not required.) This pro- 
tocol minimizes the average time one has to wait to reach 
a given purity. The decrease in this average waiting time 
over the previous protocol is a factor of two, and in this 
case the evolution of the purity is stochastic. 

In the next section we examine homodyne detection of 
a single optical qubit, and derive a deterministic rapid- 
puritifcation protocol equivalent to the first protocol dis- 
cussed above. In Section|TTT]we calculate the performance 
of two protocols that are analogous in various ways to the 
Wiseman- Ralph protocol. Section IIVI summarizes with 
some concluding remarks. 
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II. RAPID PURIFICATION FOR OPTICAL 
HOMODYNING 

The dynamics of a single mode of an optical cavity, 
where the output light is monitored via homodyne de- 
tection, is given by [13 1 

dp = -jD[a]pdt+ ^/2^{ae t0 p + pa^e^ 9 )dW 

-^(ae ie + a)er l9 )pdW, (2) 

where D[a]p = a^ap + pa^a — 2apa^ , p is the state of the 
mode, a is the mode annihilation operator, 7 is the decay 



rate of the mode from the cavity, and 77 is the efficiency of 
the photodetectors. here we have moved into the interac- 
tion picture, and thus eliminated the mode Hamiltonian 
Ho = huja)a. In this case the observer's measurement 
record is given by dr = (a + a^)dt + dW/y/Wy. If the state 
of the mode has no more than one photon, then we can 
replace a with the Pauli lowering operator <r_ = a x — ia y , 
and the SME becomes 



dp = -jD[a„]pdt + y/Zryy [<J-e ie p + pa+e ie - (a x cos6» + a y sin6>)p] dW. (3) 

We now rewrite this equation using the Bloch-sphere representation of the density matrix, a = (x, y, z), where 
p = (l/2)(/ + a • <x) and er = (a x , a y , a z ) and / is the two-by-two identity matrix. This gives 

dx = —jxdt + \J2rpi [(1 + z) cos 9 — x(xcos9 + ysin0)] dW, (4) 

dy = —jydt + 2777 [(1 + 2;) sin 6 — y{x cos 6 + y sin 9)} dW, (5) 

dz = -2j(l + z)dt- y / 2^(l + z)[x cos9 + y sin 9} dW. (6) 

Defining the "linear entropy", L, by L = 1 — Tr[p 2 ], and using the above equations we hnd that 

dL = -7{2L[1 - 77 (a; cos 9 + yarn 6)] + (77 - 1)(1 + zf) dt + ^/8rr/L (xcos9 + y sm6) dW. (7) 

I 



We wish to maximize the rate of decay of L by adjusting 
the phase of the local oscillator, 9, as the measurement 
proceeds. Inspection of the above equation makes it clear 
how to do this: we simply need to choose 9 at each time 
so that x cos 9 + ysin# = 0. This not only maximizes 
the rate of decay of L, but also eliminates the stochastic 
terms in dL and dz so that the evolutions of both are 
deterministic. This parallels the behavior of the rapid- 
purification algorithm in [l[. When we choose 9 at each 
time to maximize the rate of reduction of L, the evolution 
of L becomes 



dL 
~dt 



- 7 [2L+(7 ? -l)(l + z) 5 



(8) 



To achieve this we must continually adjust 9 so that 
9{t) = arg[y(/j) — ix(t)]. With this choice of 9 the equa- 
tion for z is simply dz/dt — —27(1 + z). We now take 
the initial state to be the maximally mixed single-qubit 
state p(0) = 1/2. Solving for the evolution of z in this 
case we have 



z(t) 



-2jt 



1. 



(9) 



and the equation of motion for the linear entropy becomes 
-7[2L+(77-l) e - 4 ^]. (10) 



dL 
~dt 



Thus the evolution of the linear entropy, under the rapid- 
purification feedback algorithm is 



-2ft 



1 1/ x / 



-2 7 n 



(11) 



We now need to compare this with the evolution of 
the average value of the linear entropy in the absence 
of any feedback. (That is, when 9 is fixed during the 
measurement.) Since we are treating the case when the 
initial state is maximally mixed, all choices for the fixed 
value of 9 are equivalent, and so we will choose 9 = for 
simplicity. When 9 is fixed the evolution of L is stochas- 
tic, and thus more complex. Nevertheless, for perfectly 
efficient detection (77 = 1) the master equation Eq.© 
is readily solved by using the linear form of the equiv- 
alent stochastic Schrodinger equation (SSE), being [l4[ 
(see also [O, [H, [3 ) , 



(12) 



d\ip)= -jcr+a-dt + ^2ja-dW \ip) 



The solution is 



P(t) 



V(t)p(0)VW 
Tr[V(t)tV(t)p(0)Y 



(13) 
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where 



V(t) 



- 7 i 



(e- 7 V + er_ + cr_a + )(l + i?cr_) 



(14) 



and R is a random variable whose probability density at 
time t is 

P(R,t) = Tr[V(tyV(t)p(0)] e ~ R J-\ (15) 



where we have defined k = (1 — e~ 27 *). When the initial 
state is the single-qubit maximally mixed state, p(0) = 
1/2, the solution is 

[e-^a+a- + (<r+ + cr_) + (l + J? 2 ) 



P(t) 
and 



i? 2 -K 



(2 + i? 2 -K) p2 /( - 9 N 



(16) 
(17) 



The evolution of the average value of the linear entropy 
is then given by 



(£(*)> 



(1 -Ti[p{t) 2 ])P{R,t)dR. (18) 



This integral cannot be solved analytically, and we will 
therefore evaluate it numerically 

When r\ is less than unity the SME is no longer equiv- 
alent to an SSE because it can increase the entropy of an 
initially pure state. Nevertheless, it turns out that it is 
possible to obtain an analytic solution to the SME by us- 
ing the above technique of solving a linear SSE. As far as 
we know this method has not appeared in the literature 
to date, and so we describe it in detail in the appendix. 

When the initial state is p(0) = 1/2 the solution for 
abitrary r] is 



P (t) = {2Afy 1 [e 



a+a- + e 



-it 



(<y+ + <y-) 



+ (1 + R 2 + [1 - rj\n) 0--0-+] , (19) 

where J\f = l+i? 2 /2 — rjn/2 is the normalization constant. 
The probability density for the random variable R is now 



P(R,t) = 



2 + R 2 - 7]K 



(20) 



We define the speed-up afforded by the rapid purifica- 
tion as the ratio of two times, s = t m /£fb- The first time, 
t m , is that taken for (L(t)) to reach a given target value 
in the absence of feedback, and the second time, ift,, is 
that taken for Lfb(t) to reach the same target value. Us- 
ing the expressions for the linear entropy in the two cases 
(Eqs. (jlip and ([18])) we plot this speed-up as a function 
of the target entropy in Figure [U and for three values of 
the detection efficiency r\. We see that in the present case 
the speed-up factor reaches a peak and then decays back 
to unity as time increases. This is quite different behav- 
ior to that of the equivalent protocol for a measurement 
of an observable (a non-dissipativc measurement), which 
tends to its maximum value as t — > oo. 



1.14 



1.12- 



1.1 



3 



1.06 



in 



1.04- 



1.02 



10 ' 



10 



10" 



Final Linear Entropy 



FIG. 1: The speed-up factor in the time required to achieve a 
given final value of the average linear entropy, (L) , afforded by 
the deterministic rapid-purification algorithm when the initial 
state of the optical qubit is completely mixed, as a function 
of (L). The various curves correspond to different values of 
the measurement efficiency r\. Solid line: 17 = 1; Dashed line: 
j] = 0.8; Dash-dot Line: r\ = 0.5. 



III. ANALOGUES OF THE WISEMAN-RALPH 
RAPID-PURIFICATION PROTOCOL 

The Wiseman-Ralph protocol minimizes the average 
time, (T), taken to reach a given linear entropy, where 
the average is once again taken over all possible mea- 
surement records 0, Q . In the case of a non-dissipative 
measurement on a single qubit, this is achieved by ap- 
plying no feedback so long as there is no Hamiltonian 
evolution. We first ask, therefore, what is the speedup 
in (T) when we apply no feedback during optical homo- 
dyning over that provided by the deterministic feedback 
algorithm derived above? In this case we resort to per- 
forming a numerical simulation of the measurement pro- 
cess (without feedback), and compare the resulting (T) 
(obtained by averaging over approximately five thousand 
trajectories) with that for the deterministic protocol as 
a function of the final linear entropy, L. For the latter 
one has a simple analytic expression for the time taken 
to reach the linear entropy L, given by solving Eq. ljlip 
for t, and this is 



T = -J-ln 

2 7 



1 - v/2 
l — i] 




(21) 



were for r\ = 1 this reduces to T = — ln(2L)/(2^). 

We plot the resulting speed-up factor in Figure for 
three values of the measurement efficiency rj. As expected 
there is a speedup, although once again this speed-up 
vanishes in the long-time limit. 

The absence of feedback is not, however, the equiva- 
lent of the Wiseman-Ralph protocol, because in our case 
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FIG. 2: The speed-up in the average time taken to reach a 
given linear entropy L for no feedback over that of the de- 
terministic feedback algorithm presented in Section [TT] The 
various curves correspond to different values of the measure- 
ment efficiency 77. Solid line: 77 = 1; Dashed line: 77 = 0.8; 
Dash-dot Line: 77 — 0.5. 



FIG. 3: The speed-up in the average time taken to reach a 
given linear entropy L for an analogue of the Wiseman- Ralph 
feedback protocol over that of the deterministic feedback al- 
gorithm presented in Section [TT] The two curves correspond 
to different values of the measurement efficiency 77. Solid line: 
77 — 1; Dashed line: 77 = 0.95. 



it does not necessarily provide the minimum value of (T). 
Without loss of generality let us choose the phase of the 
local oscillator to be 9 = 0, which is the homodyne equiv- 
alent of measuring a qubit in the a x basis. Examining the 
equations of motion for z and x in this case we see that in 
the absence of feedback the Bloch vector does not remain 
along the :c-axis, as would be true of the W-R protocol 
for a g x measurement, but in the x — z plane. This is, 
however, the closest we can come to the W-R protocol 
if we apply feedback only to the phase of the local os- 
cillator, 9; it keeps the x — y component of the Bloch 
vector aligned with the measured quadrature. Modifying 
the local oscillator is much simpler to implement experi- 
mentally that applying unitary operations to the optical 
mode, and for this reason is of most interest to us here. 

If we could implement feedback consisting of unitary 
operations on the state of the optical mode, then we could 
rotate the qubit during the measurement so as to keep the 
Bloch vector aligned with the basis of the measurement. 
For 9 — the measurement basis is the x-basis, so the 
result is z = 0, y = and the evolution is given by a 
single differential equation for x, being 



In the limit L — > (equivalcntly, t — > 00) this reduces 
to (T) m -ln(L/2)/(47). Using this expression for (T), 
and that given in Eq. (|2ip , we can immediately calculate 
the speedup factor for this protocol over the protocol in 
Section HU This is 



lim (T/(T» 

t — >OQ 



1/2. 



(24) 



Thus, unlike the previous two protocols, the speed-up 
factor for this protocol does not decay to unity as t — ► 00. 

Since we cannot obtain an analytic solution for 77 < 1, 
we calculate the speed-up factor for 77 = 0.8 by obtain- 
ing (T) using a numerical simulation and averaging over 
approximately five thousand trajectories. We plot the 
speed-up factor as a function of the final linear entropy 
for three values of 7/ in Figure [3] From this we see that 
the speed-up for this algorithm is quite a lot more sen- 
sitive to the measurement efficiency than the previous 
two. 



IV. CONCLUSION 



dx = — (1 — r})jxdt 



x 2 )dW. (22) 



When 77 = 1 the deterministic term vanishes, and the 
equation is essentially identical to that describing the 
Wiseman-Ralph protocol. In this case one can obtain 
an analytic solution for the average time to reach a given 
linear entropy L 0, [ijj • This gives 



<T) = ^H^ln 

v 1 4 7 



1 + V1-2X 



1 - VI - 2L 



(23) 



We have considered applying feedback to the homo- 
dyne detection of a single optical qubit so as to change 
the rate at which the system is purified (so called "rapid- 
purification feedback algorithms" ) . We have shown that 
there exists a feedback algorithm that increases the rate 
at which the average 'purity increases, and like its non- 
dissipative counterpart this results in a deterministic evo- 
lution for the purity of the system. Unlike its non- 
disspativc analogue, the speed-up provided by this pro- 
tocol reaches its maximum value at a finite time, decay- 
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ing to unity as t — > oo. We also found that the speed-up 
remains for measurement efficiencies well below unity, al- 
though the speed-up decreases as the efficiency drops. 

We have also examined the behavior of the average 
time taken to reach a fixed purity for different feedback 
algorithms. We found, as is true for non-dissipative mea- 
surements, that the protocol presented in Section |n] is 
slower than without feedback, as is the case for non- 
dissipative measurements. We also pointed out that the 
measurement without feedback is not the closest ana- 
logue of the Wiseman-Ralph feedback protocol for ho- 
modyne detection. We considered a more closely analo- 
gous algorithm in which the Bloch vector is aligned with 
the effective direction of the measurement, although this 
feedback cannot be implemented by merely changing the 
phase of the local oscillator. We found that this protocol 
behaves much more like those for a non-dissipative mea- 
surement, in that the speed-up factor increases monoton- 
ically and tends to a value of two in the long-time limit. 
This protocol is, however, more sensitive to noise than 
the previous protocols. The above results show that it 
should be feasible to demonstrate rapid-purification pro- 
tocols in an optical setting. 
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where dW and dV arc independent Gaussian noise 
sources so that dWdV = 0. The observer has access to 
the measurement record corresponding the measurement 
associated with dW , and thus must ultimately average 
over dV. We obtain the evolution operator which solves 
this equation by using the method given in reference [l6[ , 
and this is 



V{t,R,Q) 



e -7a f at £ Ka 2 e aR & aQ 



where 



R 



^ / e' 2 ' 1s dW{s) 
Jo 

rt 



(A2) 



(A3) 



(A4) 



and k = (1 — e~ 27 *). The probability densities for R 
and Q resulting from the above stochastic integrals are 
Gaussian, with mean zero and variances Vr — tjk and 
Vq = (1 — t])k. We will denote these Gaussian densities 
by G(R) and H{Q), respectively. 

For an initial state p(Q), the solution is thus 



Q = ^2(1-77)7 / e- 2 -> s dV(s) 



P(t, R, Q) = 



Vp(0)Vi 



(A5) 



where AT = Tr[V'Vp(Q)] is the normalization. The true 
joint probability density for R and Q is given by the prod- 
uct of the Gaussian densities G(R) and H(Q), multiplied 
by Af. That is 



APPENDIX A: SOLVING THE SME FOR 
INEFFICIENT DETECTION 

We first note that a master equation that describes in- 
efficient detection is equivalent to a master equation con- 
taining two simultaneous measurements, where the ob- 
server has access to only one, and must average over the 
results of the other [Tl| . Our method is then to solve the 
Stochastic Schrodinger equation equivalent to the SME 
with two measurements (by using the method of linear 
quantum trajectories [HI, 13, HI, [3]), and then take the 
average over the second measurement at the end to ob- 
tain the solution for inefficient detection. It turns out 
that the resulting integrals arc straightforward and give 
a fully analytic solution. The SME Eq.([2]) is thus equiv- 
alent to the linear SSE [lfj|) 



d\tp) = —ja^adt + 



y/2(l - rfjjadV 



(Al) 



P(R, Q, t) = Tr[V*Vp(0)}G(R)H(Q). (A6) 

To obtain the solution to the inefficient SME we must 
average over the Q keeping R fixed. This solution is 
therefore 

/oo 
P (R,Q,t)P(Q\R)dQ 
-OO 

1 f°° 

p(R,Q,t)P{R,Q)dQ 



P(R) 



pW)Il vmviHiQ)dQ 

1 I" 00 

— J Vp(0)V^H(Q)dQ (A7) 



where M. is merely the normalization. From this we see 
that we need only perform an integration over the Gaus- 
sian density for Q, which is straightforward. 
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